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Vacuum polarization integrals involve the vector spectral functions which can be experimentally determined 
from two sources: (i) e'^e~ annihilation cross sections and (ii) hadronic r decays. Recently results with compara- 
ble precision have become available from CMD-2 on one side, and ALEPH, CLEO and OPAL on the other. The 
comparison of the respective spectral functions involves a correction from isospin-breaking effects, which is eval- 
uated. After the correction it is found that the dominant tttv spectral functions do not agree within experimental 
and theoretical uncertainties. Some disagreement is also found for the Att spectral functions. The consequences 
of these discrepancies for vacuum polarization calculations are presented, with the emphasis on the muon anoma- 
lous magnetic moment. Adding quadratically experimental and theoretical uncertainties, we find that the full 
Standard Model prediction of deviates from the recent BNL measurement at the level of 3.0 (e'''e~-based) and 
0.9 (r-based) standard deviations. 



1. INTRODUCTION 

Hadronic vacuum polarization in the photon 
propagator plays an important role in the preci- 
sion tests of the Standard Model. This is the case 
for the evaluation of the electromagnetic coupling 
at the Z mass scale, q;(M^), which receives a con- 
tribution AQ;had(A'/|) of the order of 2.8 10"^ 
that must be known to an accuracy of better 
than 1% so that it does not limit the accuracy 
on the indirect determination of the Higgs boson 
mass from the measurement of sin^ 9\y . Another 
example is provided by the anomalous magnetic 
moment = (g^ — 2)/2 of the muon, where the 
hadronic vacuum polarization component is the 
leading contributor to the uncertainty of the the- 
oretical prediction. 

Starting from Refs. [1,2] there is a long history 
of calculating the contributions from hadronic 
vacuum polarization in these processes. As they 
cannot be obtained from first principles because 
of the low energy scale involved, the computation 
relies on analyticity and unitarity so that the rele- 



vant integrals can be expressed in terms of an ex- 
perimentally determined spectral function which 
is proportional to the cross section for e~^e~ anni- 
hilation into hadrons. The accuracy of the calcu- 
lations has therefore followed the progress in the 
quality of the corresponding data [3]. Because 
the latter was not always suitable, it was deemed 
necessary to resort to other sources of informa- 
tion. One such possibility was the use [4] of the 
vector spectral functions derived from the study 
of hadronic r decays [5] for the energy range less 
than 1.8 GeV. Another one occurred when it was 
realized in the study of r decays [6] that pertur- 
bative QCD could be applied to energy scales as 
low as 1-2 GeV, thus offering a way to replace 
poor e+e" data in some energy regions by a reli- 
able and precise theoretical prescription [7-9] . Fi- 
nally, without any further theoretical assumption, 
it was proposed to use QCD sum rules [10,11] in 
order to improve the evaluation in energy regions 
dominated by resonances where one has to rely 
on experimental data. Using these improvements 
the lowest-order hadronic contribution to was 
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found to be [11] 

^had.LO ^ (692.4 ±6.2) 10^1° 



(1) 



The complete theoretical prediction includes in 
addition QED, weak and higher order hadronic 
contributions. 

The anomalous magnetic moment of the muon 
is exporinientally known to very high accuracy. 
Combined with the older, less precise results from 
CERN [12], the measurements from the E821 ex- 
periment at BNL [13-15], including the most re- 
cent result [16], yield 

a^^P = (11 659 203 ±8) 10"^° , (2) 

and are aiming at an ultimate precision of 4 10~^° 
in the future. The previous experimental re- 
sult [15] was found to deviate from the theo- 
retical prediction by 2.6 a, but a large part of 
the discrepancy was originating from a sign mis- 
take in the calculation of the small contribution 
from the so-called light-by-light (LBL) scattering 
diagrams [17,18]. The new calculations of the 
LBL contribution [19-21] have reduced the dis- 
crepancy to a nonsignificant 1.6 cr level. At any 
rate it is clear that the presently achieved exper- 
imental accuracy already calls for a more precise 
evaluation of a^.'^'^'-'. 

New experimental and theoretical develop- 
ments have prompted the re-evaluation of the 
hadronic contributions presented in Ref. [22] and 
reported here: 

• new, precise results have been obtained at 
Novosibirsk with the CMD-2 detector in the 
region dominated by the p resonance [23], 
and more accurate R measurements have 
been performed in Beijing with the BES de- 
tector in the 2-5 GeV energy range [24]. 

• new preliminary results are available from 

the final analysis of r decays with ALEPH 
using the full statistics accumulated at 
LEPl [25]; also the information from the 
spectral fimctions measured by CLEO [26, 
27] and OPAL [28] has been incorporated 
in the analysis. 

• new results on the evaluation of isospin 
breaking have been produced [29-31], thus 



providing a better understanding of this 
critical area when relating vector r and 
isovector e+e~ spectral functions. 

2. MUON MAGNETIC ANOMALY 

It is convenient to separate the Standard Model 
prediction for the anomalous magnetic moment of 
the muon into its different contributions, 



„SM _ „QED I had , weak 



(3) 



with aji^'i = 0,^'^'^° + aji^'i'^o + a^^'^^^ , 
where a^^^ = (11658 470.6 ± 0.3) 10"!° is the 
pure electromagnetic contribution (see [32,33] 
and references therein), a^ad.LO jg ^j^g jowest- 
order contribution from hadronic vacuum polar- 
ization, a^"'^'^'^ = (-10.0±0.6) 10-1" is the cor- 
responding higher-order part [34,4], and aj^^^^ = 
(15.4 ± 0.1 ± 0.2) 10-1°, where the first error 
is the hadronic uncertainty and the second is 
due to the Higgs mass range, accounts for cor- 
rections due to exchange of the weakly interact- 
ing bosons up to two loops [35]. For the LBL 
part we add the values for the pion-pole contribu- 
tion [19-21] and the other terms [20,21] to obtain 

^had,LBL ^ (8 g ^ 3 iQ-W 

By virtue of the analyticity of the vacuum 
polarization correlator, the contribution of the 
hadronic vacuum polarization to can be cal- 
culated via the dispersion integral [36] 



a 



had.LO 



Q^(0) 
3^2 



ds ^R{s) , 



(4) 



4mJ 



where K{s) is the QED kernel [37] strongly em- 
phasizing the low-energy spectral functions in the 
integral (4). In effect, about 91% of the total 



contribution to 



had.LO 



is accumulated at center- 



of-mass encirgies \/s below 1.8 GeV and 73% of 
i,^had.LO (;overed by the two-pion final state 
which is dominated by the p(770) resonance. In 
Eq. (4), R{s) = i?(°)(s) denotes the ratio of the 
'bare' cross section for e+e- annihilation into 
hadrons to the pointlike muon-pair cross section. 
At low mass-squared, R{s) is taken from experi- 
ment. 
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3. DATA FROM e+e" ANNIHILATION 

The exclusive low energy e~^e~ cross sections 
have been measured mainly by experiments run- 
ning at e^e~ colliders in Novosibirsk and Orsay. 
Due to the high hadron multiplicity at energies 
above ~ 2.5 GeV, the exclusive measurement of 
the respective hadronic final states is not practi- 
cable. Consequently, the (experiments at the high 
energy colliders have measured the total inclusive 
cross section ratio R. 

The most precise data from CMD-2 on the 
e'^e~—^ 7r+7r~ cross sections are now available in 
their final form [23] . They differ from the prelim- 
inary ones, released two years ago [38], mostly in 
the treatment of the radiative corrections. The 
various changes resulted in a reduction of the 
cross section by about 1% below the p peak and 
5% above. The overall systematic error of the fi- 
nal data is quoted to be 0.6% and is dominated 
by the uncertainties in the radiative corrections 
(0.4%). Agreement is observed between CMD- 
2 and the previous experiments within the much 
larger uncertainties (2-10%) quoted by the latter. 

Large discrepancies are observed between the 
different data sets for the e~^e~—^ tt+tt" 77*^77*^ 
cross sections (see Fig. 5). These are probably 
related to problems in the calculation of the de- 
tection efficiency, since the efficiencies are small 
in general (~ 10 — 30%) and they are affected by 
uncertainties in the decay dynamics that is as- 
sumed in the Monte Carlo simulation. One could 
expect the more recent experiments (CMD-2 [40] 
and SND [41]) to be more reliable in this context 
because of specific studies performed in order to 
identify the major decay processes involved. 

For e+e~ ^ 7r+7r~7r+7r~ the experiments agree 
reasonably well within their quoted uncertainties 
(see Fig. 4 in Section 7). 

A detailed compilation and complete references 
of all the data used for this analysis are given in 
Ref. [22]. 



ALEPH [25,5], CLEG [26,27] and OPAL [28]. 
The branching fraction B^j^o for the r — > 
Vr 7r~7r° (7) decay mode is of particular interest 
since it provides the normalization of the corre- 
sponding spectral function. The new value [25], 
B^^o = (25.47 ± 0.13) %, turns out to be larger 
than the previously published one [42] based on 
the 1991-93 LEPl statistics, (25.30 ± 0.20) %. 

In the limit of isospin invariance, the corre- 
sponding e+e~ isovector cross sections are cal- 
culated via the isospin rotations 



(5) 



„ 47ra^ 

2 v^- 3^0 , (6) 



47ra^ 



^^27r-7r+7rO 



■ SttO 



(7) 



(8) 



The r spectral function vv{s) for a given vector 
hadronic state V is defined by [43] 



vv{s) 



(X 



Vt e~ Ue) Nv ds 



B{t- 



mt 



(9) 



where = 0.9748 ± 0.0010 (using [44]). The 

spectral functions are obtained from the cor- 
responding invariant mass distributions, after 
subtracting out the non-T background and the 
feedthrough from other r decay channels, and af- 
ter a final unfolding from detector effects such as 
energy and angular resolutions, acceptance, cal- 
ibration and photon identification. Agreement 
within errors is observed for the three input spec- 
tral functions v^-^o (ALEPH, CLEG, GPAL), so 
that we use their weighted average in the follow- 
ing. 



4. DATA FROM HADRONIC r DECAYS 

Data from t decays into two- and four- 
pion final states r~ — > i/r7r~7r°, t~ — > Z/'T-7r~37r° 
and T~ — > i/r27r~7r+7r'^, are available from 



5. RADIATIVE CORRECTIONS FOR 
e+e- DATA 

Radiative corrections applied to the measured 
e+e~ cross sections are an important step in the 
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experimental analyses. They involve the consid- 
eration of several physical processes and lead to 
large corrections. We stress that the evaluation 
of the integral in Eq. (4) requires the use of the 
'bare' hadronic cross section, so that the input 
data must be analyzed with care in this respect. 
Several steps are to be considered: 

• Corrections are applied to the luminos- 
ity determination, based on large-angle 
Bhabha scattering and muon-pair produc- 
tion in the low-energy experiments, and 
small-angle Bhabha scattering at high en- 
ergies. These processes are usually cor- 
rected for external radiation, vertex correc- 
tions and vacuum polarization from lepton 
loops. 

• The hadronic cross sections given by the ex- 
periments are in general corrected for initial 
state radiation and the effect of loops at the 
electron vertex. 

• The vacuum polarization correction in 

the photon propagator is a more delicate 
point. The cross sections need to be 
fully corrected for our use, i.e. fTbare = 
a'drosscd(a(0)/a(s))2, where (Tdrcsscd is the 
measured cross section already corrected for 
initial state radiation, and a(s) is obtained 
from rcsummation of the lowest-order eval- 
uation, giving a{s) = a{0)/{l — Aaicp{s) — 
Aahad(s))- Whereas AQ;iep(s) can be cal- 
culated analytically. Aahad(s) is related by 
analyticity and unitarity to a dispersion in- 
tegral, akin to Eq. (4). Since the hadronic 
correction involves the knowledge of R{s) at 
all energies, including those where the mea- 
surements are made, the procedure has to 
be iterative, and requires experimental as 
well as theoretical information over a large 
energy range. 

The new data from CMD-2 [23] are 
explicitly corrected for both leptonic 
and hadronic vacuum polarization effects, 
whereas the preliminary data from the same 
experiment [38] and data from other exper- 
iments were not (see follow up of this dis- 
cussion in Ref. [22]) 



Table 1 

Expected sources of isospin symmetry breaking 

between e~^e~ and t spectral functions in the 27r 
and 47r channels, and the corresponding correc- 
tions to a}}.'^'^'^ as obtained from r data. 



Sources of Isospin Aa^^'^^ (10-^°) 

Symmetry Breaking tt+tt" 



Short distance rad. corr. —12.1 ± 0.3 

Long distance rad. corr. —1.0 

m-TT- 7^ TO^o (/3 in cross section) —7.0 

m^- ^ m-j^o {(3 in p width) -1-4.2 

nip- ^ nipo ± 2.0 

p — oj interference -1-3.5 ± 0.6 

Electromagnetic decay modes —1.4 ± 1.2 



Sum -13.8 ±2.4 



• In Eq. (4) one must incorporate in R{s) 
the contributions of all hadronic states pro- 
duced at the energy ^/s. In particular, ra- 
diative effects in the hadronic final state 
must be considered, i.e., final states such 
asV + j (FSR) have to be included. While 
FSR has been added to the newest CMD- 
2 data [23], this is not the case for the 
older data and thus has to be corrected 
for [22] "by hand", using an analytical ex- 
pression computed in scalar QED (point- 
like pions) [45]. 

In summary, we correct each e+e~ experimen- 
tal result, but those from CMD-2 (tttt), by the 
factor Chvp • C'fsR) where Chvp accounts for 
hadronic vaccum polarization and CpsR stands 
for the FSR correction. We assign uncertainties of 
50% (vacuum polarization) and 100% (FSR cor- 
rections), which are considered to be fully corre- 
lated between all channels to which the correc- 
tions apply. 

6. ISOSPIN BREAKING IN e+e" AND r 
SPECTRAL FUNCTIONS 

The relationships (5), (6) and (7) between e+e" 
and T spectral functions only hold in the limit of 
exact isospin invariance. This is the Conserved 
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Vector Current (CVC) property of weak decays. 
It follows from the factorization of strong inter- 
action physics as produced through the 7 and W 
propagators out of the QCD vacuum. However, 
we know that we must expect symmetry break- 
ing at some level from electromagnetic effects and 
even in QCD because of the up and down quark 
mass splitting. Since the normalization of the r 
spectral functions is experimentally known at the 
0.5% level, it is clear that isospin-breaking effects 
must be carefully examined if one wants this pre- 
cision to be maintained in the vacuum polariza- 
tion integrals. 

Because of the dominance of the tttt contribu- 
tion in the energy range of interest for r data, 
we discuss mainly this channel, following our ear- 
lier analysis [4]. The corrections on aj^*^'^*-* from 
isospin breaking are given in Table 1. A more 
complete discussion, in particular with respect 
to the corrections previously applied is given in 
Rcf. [22]. 

The dominant contribution to the electroweak 
radiative corrections stems from the short dis- 
tanc;c correction to the effective four-fcrmion cou- 
pling T~ Vr{du)~ enhancing the r amplitude 
by the factor S^^ = 1.0194 [46]. This correction 
leaves out the possibility of sizeable contributions 
from virtual loops. This problem was studied in 
Ref. [30] within a model based on Chiral Pertur- 
bation Theory. In this way the correct low-energy 
hadronic structure is implemented and a consis- 
tent framework has been set up to calculate elec- 
troweak and strong processes, such as the radia- 
tive corrections in the r — » VtTT~-k^ decay. Their 
new analysis [31] directly applies to the inclusive 
radiative rate, r i/T-7r~7r*' (7), as measured by 
the experiments. The relation between the Born 
level e+e~ spectral function and the r spectral 
function reads [31] 

2 



1 PI 

Gem 



where Gem = Gem(s) is the long-distance radia- 
tive correction involving both real photon emis- 
sion and virtual loops (the infrared divergence 
cancels in the sum). Note that the short-distance 
Sew correction, discussed above, is already ap- 
plied in the definition of V- {cf. Eq. (9)), but 



its value differs from the one given in Ref. [46], 
because subleading quark-level and hadron-level 
contributions should not be added, as double 
counting would occur. The correct expression for 
the 7r~7r° mode therefore reads 



'-'ew'-'EM 

Qsub,lep 



(1.0233 ±0.0006) -Gem 



(11) 



the subleading hadronic corrections being now in- 
corporated in the mass-squarcd-depcndent Gem 
factor. Equation (10) explicitly corrects for the 
mass difference between neutral and charged pi- 
ons affecting the cross section and the width of 
the p. 

The different contributions to the isospin- 
breaking corrections are shown in the second col- 
umn of Table 1. The dominant uncertainty stems 
from the p^-p^ mass difference. 

Since the integral (4) requires as input the 
e+e~ spectral function including FSR photon 
emission, a final correction is necessary. It is iden- 
tical to that applied in the CMD-2 analysis [23,45] 
{cf. Section 5). The total correction to the t re- 
sult amounts to (-9.3 ± 2.4) 10"^°. 

There exists no comparable study of isospin 
breaking in the 47r channels. Only kinematic 
corrections resulting from the pion mass differ- 
ence have been considered so far [29], which we 
have applied in this analysis. It creates shifts of 
-0.7 10-1° (-3.8%) and +0.1 10-^° {+1.1%) for 
27r"'"27r- and ir'^ ■jt~ 2n'^ , respectively. 

7. COMPARISON OF e+e" AND r 
SPECTRAL FUNCTIONS 

The e+g- and the isospin-breaking corrected r 
spectral functions can be directly compared for 
the dominant 27r and 47r final states. For the 27r 
channel, the p-dominated form factor falls off very 
rapidly at high energy so that the comparison 
can be performed in practice over the full energy 
range of interest. The situation is different for the 
47r channels where the r decay kinematics limits 
the exercise to energies less than ~ 1.6 GeV, with 
only limited statistics beyond. 

Fig. 3 shows the comparison for the 27r spectral 
functions. Visually, the agreement seems satisfac- 
tory, however the large dynamical range involved 
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Figure 3. Comparison of the tt+tt^ spectral functions from e^e^ and isospin-breaking corrected r data, 
expressed as e^e~ cross sections. The band indicates the combined e^e~ and r result within Icr errors. 
It is given for illustration purpose only. 



7 



T 




0.25 0.5 0.75 

s (GeV^) 



Figure 1. Relative comparison of the tt+tt^ spec- 
tral functions from e"'"e~ and isospin-breaking 
corrected t data, expressed as a ratio to the r 
spectral function. The band shows the uncer- 
tainty on the latter. 



T 



. • CMD-2 

T Average 




0.5 0.55 0.6 0.65 0.7 

s (GeV^) 

Figure 2. Relative comparison in the p region 
of the 7r"'"7r~ spectral functions from e^e~ and 
isospin-breaking corrected r data, expressed as a 
ratio to the r spectral function. The band shows 
the uncertainty on the latter. 



does not permit an accurate test. To do so, the 
e'^e^ data are plotted as a point-by-point ratio 
to the T spectral function in Fig. 1, and enlarged 
in Fig. 2, to better emphasize the region of the p 
peak. The e+e^ data are significantly lower by 
2-3% below the peak, the discrepancy increasing 
to about 10% in the 0.9-1.0 GeV region. 

The comparison for the 47r cross sections is 
given in Fig. 4 for the 27r+27r~ channel and in 
Fig. 5 for 7r+7r^27r'^. The latter suffers from large 
differences between the results from the different 
e^e~ experiments. The r data, combining two 
measured spectral functions according to Eq. (7) 
and corrected for isospin breaking as discussed in 
Section 6, lie somewhat in between with large un- 
certainties above 1.4 GeV because of the lack of 
statistics and a large feedthrough background in 
the T ^ Vr 7r~37r° mode. In spite of these difficul- 
ties the 7r~37r° spectral function is in agreement 
with e+e~ data as can be seen in Fig. 4. It is 
clear that intrinsic discrepancies exist among the 
e^e~ experiments and that a quantitative test of 
CVC in the 7r+7r^27r*' channel is premature. 



7.1. Branching Ratios in r Decays and 
CVC 

A convenient way to assess the compatibil- 
ity between e^e~ and r spectral functions pro- 
ceeds with the evaluation of t decay fractions us- 
ing the relevant e^e~ spectral functions as in- 
put. All the isospin-breaking corrections dis- 
cussed in Section 6 are included. The advan- 
tage of this procedure is to allow a quantitative 
comparison using a single number. The weight- 
ing of the spectral function is however differ- 
ent from the vacuum polarization kernels. Us- 
ing the branching fraction B{t^ — > Vr e.~ Ve) = 
(17.810 ± 0.039)%, obtained assuming leptonic 
universality in the charged weak current [25] , the 
results for the main channels are given in Ta- 
ble 2. The errors quoted for the CVC values 
are split into uncertainties from {i) the experi- 
mental input (the e+e~ annihilation cross sec- 
tions) and the numerical integration procedure, 
{ii) the missing radiative corrections applied to 
the relevant e^e^ data, and {iii) the isospin- 
breaking corrections when relating r and e^e~ 
spectral functions. The values for the r branch- 
ing ratios involve measurements [25,47,48] given 
without charged hadron identification, i.e., for 
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Table 2 

Branching fractions of r vector decays into 2 and 4 pions in the final state. Second column: world 

average. Third column: inferred from e+e~ spectral fimctions using the isospin relations (5-7) and 
correcting for isospin breaking. The experimental error of the tt+tt" CVC value contains an absolute 
procedural integration error of 0.08%. Experimental errors, including uncertainties on the integration 

procedure, and theoretical (missing radiative corrections for e+e~, and isospin-brcaking corrections and 
Vud for t) arc shown s(^parat(^Iy. Right coliunn: differences Ik^avch^u the? direct nicasurc^mcHits in r decays 
and I1k> CVC (valuations, where the seijurate terrors liaA'c^ Ixnni ackknl in ([iiudrat ur(\ 

■|^^^^ Branching fractions (in %) 

°^ T data e+e" via CVC A(T-e+e") 

r-^z/^7r-7r° 25.46 ±0.12 23.98 ± 0.25exp ± O.llrad ± 0.12su(2) +1.48 ±0.32 



0.30 

r- ^ z^xTt-Stt" 1.01 ±0.08 1.09 ± 0.06exp ± 0.02rad ± 0.05su(2) -0.08 ±0.11 



0.08 

T- ^ Ur2'K-'K+7T° 4.54 ±0.13 3.63 ± 0.19exp ± 0.04rad ± 0.09su(2) +0.91 ±0.25 

^ V ' 

0.21 



the hiT^VTi hiiT^Vr and 'ihTr'-^iyr final states. The 
corresponding channels with charged kaons have 
been measured [49,50] and their contributions can 
be subtracted out in order to obtain the pure pi- 
onic modes. As expected from the preceding dis- 
cussion, a large discrepancy is observed for the 
T ^ 7r^7r° branching ratio, with a difference of 

(-1.48 ± 0.12^ ± 0.25oo ± O.llrad ± 0.12sXT(2))%, 

where the uncertainties are from the r branching 
ratio, e+e^ cross sections, e+e^ missing radia- 
tive corrections and isospin-brcaking corrections 
(including the uncertainty on Vud), respectively. 
Adding all errors in quadrature, the effect repre- 
sents a 4.6 a discrepancy. Since the disagreement 
between e+e~ and r spectral functions is more 
pronounced at energies above 750 MeV, we ex- 
pect a smaller discrepancy in the calculation of 
^had.LO because of the steeply falling kernel K{s) 
in this case. 

The situation in the An channels is different. 
Agreement is observed for the 7r~37r'' mode within 
an accuracy of 11%, however the comparison is 
not satisfactory for the 27r~7r+7r° mode. In the 
latter case, the relative difference is very large, 
(22 ± 6)%, compared to any reasonable level of 
isospin symmetry breaking. As such, it rather 
points to experimental problems that have to be 



investigated. 

8. SPECIFIC CONTRIBUTIONS 

In some energy regions where data information 
is scarce and reliable theoretical predictions are 
available, we use analytical contributions to ex- 
tend the experimental integral. Also, the treat- 
ment of narrow resonances involves a specific pro- 
cedure. 

8.1. The TT+TT" Threshold Region 

To overcome the lack of precise data at thresh- 
old energies and to benefit from the analyticity 
property of the pion form factor, a third order 
expansion in s is used. The pion form factor 
is connected with the tt+tt" cross section via 
— (3s/7ra^/3Q) a^TT- The expansion for small 
s reads 

= 1 + ^{r\ s + ci + C2 + O(s^) , (12) 

where we use (r^)^ = (0.439 ± 0.008) fm^ [51], 
and the two parameters ci,2 are fitted to the data 
in the range [2171^^, 0.6 GeV]. The results of the 
fits are explicitly quoted in Ref. [22]. We show 
the functions obtained in Fig. 6. Good agree- 
ment is observed in the low energy region where 
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1 1.5 2 2.5 3 3.5 4 4.5 5 



s (GeV^) 

Figure 4. Comparison of the 27r+27r~ spectral 
functions from e'^e~ and isospin-breaking cor- 
rected r data, expressed as e+e^ cross sections. 




1 1.5 2 2.5 3 3.5 



s (GeV^) 

Figure 5. Comparison of the 7r+7r~27r° spectral 
functions from e+e~ and isospin-breaking cor- 
rected r data, expressed as e+e~ cross sections. 



the expansion should be reliable. Since the fits 
incorporate unquestionable constraints from first 
principles, we have chosen to use this parameter- 
ization for evaluating the integrals in the range 
up to 0.5 GeV. 

8.2. Integration over the uj and Reso- 
nances 

In the regions around the lo and resonances 
we have assumed in the preceding works that the 
cross section of the 7r"'"7r~7r'^ production on the 
one hand, and the Tr+Tr^Tr", K^K^ as well as 
KgK^ production on the other hand is saturated 
by the corresponding resonance production. In a 
data driven approach it is however more careful to 
directly integrate the measurement points with- 
out introducing prior assumptions on the underly- 
ing process dynamics [52]. Possible non-resonant 
contributions and interference effects are thus ac- 
counted for. 

Notwithstanding, a straightforward trapezoidal 
integration buries the danger of a bias: with in- 
sufficient scan density, the linear interpolation of 
the measurements leads to a significant overesti- 



mation of the integral when dealing with strongly 
concave functions such as the tails of Breit- 
Wigner resonance curves. 

We therefore perform a phenomenological fit 
of a BW resonance plus two Gaussians (only 
one Gaussian is necessary for the uj) to account 
for contributions other than the single resonance. 
Both fits result in satisfactory values. We have 
accounted for the systematics due to the arbi- 
trariness in the choice of the parametrization by 
varying the functions and parameters used. The 
resulting effects are numerically small compared 
to the experimental errors. 

Since the experiments quote the cross sec- 
tion results without correcting for leptonic and 
hadronic vacuum polarization in the photon prop- 
agator, we perform the correction here. The cor- 
rection of hadronic vacuum polarization being it- 
erative and thus only approximative, we assign 
half of the total vacuum polarization correction 
as generous systematic errors (c/. Section 5). In 
spite of that, the evaluation of ajj'*'^'^'-* is dom- 
inated by the experimental uncertainties. Since 
the trapezoidal rule is biased, we choose the re- 
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0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 



s (GeV^) 

Figure 6. Fit of the pion form factor from 4m^ 
to 0.35 GeV^ using a third-order Taylor expan- 
sion with the constraints at s = and the mea- 
sured pion r.m.s. charge radius from space-Hke 
data [51]. The result of the fit is integrated only 
up to 0.25 GeV^. 



suits based on the BW fits for the final evaluation 
of a{;'^<^'LO. 

8.3. Narrow cc and bb Resonances 

The contributions from the narrow J/tp reso- 
nances are computed using a relativistic Breit- 
Wigner parametrization for their line shape. The 
physical values for the resonance parameters and 
their errors are taken from the latest compila- 
tion in Ref . [44] . Vacuum polarization effects are 
already included in the quoted leptonic widths. 
The total parametrization errors are then calcu- 
lated by Gaussian error propagation. This inte- 
gration procedure is not followed for the ip{3S) 
state which is already included in the R measure- 
ments, and for the T resonances which are repre- 
sented in an average sense (global quark-hadron 
duality) by the bb QCD contribution, discussed 
next. 

8.4. QCD Prediction at High Energy 

Since the emphasis in this paper is on a com- 
plete and critical evaluation of spectral functions 
from low-energy data, we have adopted the con- 
servative choice of using the QCD prediction only 
above an energy of 5 GeV. The details of the 



calculation can be found in our earlier publica- 
tions [7,11,22] and in the references therein. 

A test of the QCD prediction can be performed 
in the energy range between 1.8 and 3.7 GeV. The 
contribution to ah^-d.LO ^^^is region is computed 
to be (33.87±0.46) 10"^° using QCD, to be com- 
pared with the result, (34.9 ± 1.8) 10"^" from the 
data. The two values agree within the 5% accu- 
racy of the measurements. 

9. RESULTS 

9.1. Lowest Order Hadronic Contributions 

We use the trapezoidal rule to integrate the 
experimental data points (with the exception of 
the narrow resonances). Correlations between the 
measurements as well as among experiments have 
been taken into account [22]. Before adding up all 
the contributions to a^^'^'^'-', we shall summarize 
the procedure. On the one hand, the e^e~-based 
evaluation is done in three pieces: the sum of 
exclusive channels below 2 GeV, the R measure- 
ments in the 2-5 GeV range and the QCD pre- 
diction for R above. Major contributions stem 
from the 27r (73%) and the two Att (4.5%) chan- 
nels. On the other hand, in the r-based evalua- 
tion, the latter three contributions are taken from 
T data up to 1.6 GeV and complemented by e^e~ 
data above, because the r spectral functions run 
out of precision near the kinematic limit of the r 
mass. Thus, for nearly 77% of ah3,d,LO (contribut- 
ing 80% of the total error-squared), two indepen- 
dent evaluations {e~^e~ and t) are produced, the 
remainder being computed from e+e^ data and 
QCD alone. 

Fig. 7 gives a panoramic view of the e~^e~ data 
in the relevant energy range. The shaded band 
below 2 GeV represents the sum of the exclu- 
sive channels considered in the analysis. It turns 
out to be smaller than our previous estimate [4], 
essentially because more complete data sets are 
used and new information on the dynamics could 
be incorporated in the isospin constraints for the 
missing channels. The QCD prediction is indi- 
cated by the cross-hatched band. It is used in 
this analysis only for energies above 5 GeV. Note 
that the QCD band is plotted taking into account 
the thresholds for open flavour B states, in order 
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Figure 7. Compilation of the data contributing to aJJ ' . Shown is the total hadronic over muonic cross 
section ratio R. The shaded band below 2 GcV represents the sum of the exclusive channels considered 
in this analysis, with the exception of the contributions from the narrow resonances which are given as 
dashed lines. All data points shown correspond to inclusive measurements. The cross-hatched band gives 
the prediction from (essentially) perturbative QCD, which is found to be in good agreement with the 
measurements in the continuum above 2 GeV. In this figure the bb threshold is indicated at the onset of 
BB states in order to facilitate the comparison with data in the continuum. In the actual calculation the 
threshold is taken at twice the pole mass of the b quark. 
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Table 3 

Summary of the ajj*^'^*-* contributions from e+e~ annihilation and r decays. The uncertainties on the 
vacuum polarization and FSR corrections are given as second errors in the individual e+e^ contribu- 
tions, while those from isospin breaking are similarly given for the t contributions. These 'theoretical' 
uncertainties are correlated among all channels, except in the case of isospin breaking which shows little 
correlation between the 2tt and Att channels. The errors given for the sums in the last line are from the 
experiment, the missing radiative corrections in e+e" and, in addition for r, SU(2) breaking. 

„had,LO (-\r\-10\ 

Modes Energy [GeV] , " ^ ' , 

^ ^ e+e~ t(^) A(e+e--T) 



Low S exp. TT+TT" 


[2rn^± - 


0.500] 


58.04 ± 1.70 ± 1.14 


56.03 ± 1.61 ±0.28 


+2.0 ±2.6 


TT+TT" 


[0.500 - 


1.800] 


440.81 ± 4.65 ± 1.54 


464.03 ±3.19 ±2.34 


-23.2 ±6.3 


7r°7, 777 


[0.500 - 


1.800] 


0.93 ±0.15 ±0.01 


- 


- 


w 


[0.300 - 


0.810] 


36.94 ±0.84 ±0.80 


- 


- 


TT+TT"??" [below (t>] 


[0.810- 


1.000] 


4.20 ±0.40 ±0.05 


- 


- 


(P 


[1.000- 


1.055] 


34.80 ±0.92 ±0.64 






TT'^TT~TT^ [above </)] 


[1.055- 


1.800] 


2.45 ±0.26 ±0.03 






7r+7r-27r° 


[1.020- 


1.800] 


16.73 ± 1.32 ±0.20 


21.44 ± 1.33 ±0.60 


-4.7 ± 1.8 


27r+27r- 


[0.800 - 


1.800] 


13.95 ±0.90 ±0.23 


12.34 ±0.96 ±0.40 


+1.6 ±2.0 


27r+27r-7r° 


[1.019- 


1.800] 


2.09 ±0.43 ±0.04 






7r+7r-37rO(2) 


[1.019- 


1.800] 


1.29 ±0.22 ±0.02 






Stt+Stt" 


[1.350- 


1.800] 


0.10 ±0.10 ±0.00 






27r+27r-27r° 


[1.350- 


1.800] 


1.41 ±0.30 ±0.03 








[1.350- 


1.800] 


0.06 ±0.06 ±0.00 






Tji-^ 7r+7r~7, 27)77+77- 


[1.075- 


1.800] 


0.54 ±0.07 ±0.01 






uj{-^ 7707)77° 


[0.975 - 


1.800] 


0.63 ±0.10 ±0.01 






7r07)(7777)° 


[1.340- 


1.800] 


0.08 ±0.01 ±0.00 






K+K- 


[1.055 - 


1.800] 


4.63 ±0.40 ±0.06 








[1.097- 


1.800] 


0.94 ±0.10 ±0.01 






K°K^n^ (2) 


[1.340- 


1.800] 


1.84 ±0.24 ±0.02 






KKn° (2) 


[1.440- 


1.800] 


0.60 ±0.20 ±0.01 






KKtttt (2) 


[1.441 - 


1.800] 


2.22 ± 1.02 ±0.03 






R^J2 6xcl. modes 


[1.800- 


2.000] 


8.20 ±0.66 ±0.10 






R [Data] 


[2.000- 


3.700] 


26.70 ±1.70 ±0.00 








[3.088 - 


3.106] 


5.94 ±0.35 ±0.03 






^(25) 


[3.658 - 


3.714] 


1.50 ±0.14 ±0.00 






R [Data] 


[3.700 - 


5.000] 


7.22 ±0.28 ±0.00 






Rudsc [QCD] 


[5.000 - 


9.300] 


6.87 ±0.10 ±0.00 






Rudscb [QCD] 


[9.300 - 


12.00] 


1.21 ±0.05 ±0.00 






-R«dsc5t [QCD] 


[12.0- 


- 00] 


1.80 ±0.01 ±0.00 







E(e+e-^hadrons) [2m.± - oo] 684.7±6.0e.p 709.0±5.1 _24.3±7.9,ot 

± 3. brad ± 1-^rad ± ^•OSU(2) 

^ Not including lu and <f) resonances (see text). 

Using isospin relations (see text). 
^ e+e- data are used above 1.6 GeV (see text). 
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Figure 8. Comparison of the results with the BNL 
measurement [16]. Also shown are our previous 
estimates [3,11] obtained before the CMD-2 data 
were available, and the recent evaluation of Hagi- 
wara et al. [52]. 



to facilitate the comparison with the data in the 
continuum. However, for the evaluation of the 
integral, the bb threshold is taken at twice the 
pole mass of the b quark, so that the contribution 
includes the narrow T resonances, according to 
global quark-hadron duality. 

The contributions from the different processes 
in their indicated energy ranges are listed in 
Table 3. Wherever relevant, the two e+e^- 
and T-based evaluations are given. The dis- 
crepancies discussed above are now expressed di- 
rectly in terms of a^'^'^^ giving smaller esti- 
mates for e+e" data by (— 21.2±6.4oxp±2.4i.ad± 
2-6su(2) (±7.3totai)) 10"-^° for the 27r channel and 

(-3.1±2.6exp±0.3,ad±1.0sU(2) (±2.9total)) IQ-^" 

for the sum of the 47r channels. The to- 
tal discrepancy (—24.3 ± 6.9oxp ± 2.7rad ± 
2-8su(2) (±7.9totai)) 10^^" amounts to 3.1 stan- 
dard deviations and precludes from performing 
a straightforward combination of the two evalua- 
tions. 



9.2. Results for 

The results for the lowest order hadronic con- 
tribution are 



had.LO 
.,had,LO 



(684.7 ±6.0exp±3.6,ad) 10" 
(709.0 ±5.1cxp±3.0,a 



-10 



10 



-10 



where the first numbers are e^e^-based and the 
second r-based. Adding the QED, higher-order 
hadronic, light-by-light scattering and weak con- 
tributions as given in Section 2, we obtain for 

af^ = (11 659 169.3 ±7.0had±3.5LBL 

± 0.4qed+ew) 10-" , 

af = (11 659 193.6 ±5.9had±3.5LBL 

± 0.4qed+ew) 10"" , 

where again the first results are e"'"e~-based and 
the second r-based. These values can be com- 
pared to the present experimental average given 
in Eq. (2). Adding experimental and theoreti- 
cal errors in quadrature, the differences between 
measured and computed values are found to 
be (first number e"'"e~-based and the second t- 
based) 



tCxp 



.,SM 



= (33. 7± 11.2) 10"" 
= (9.4 ± 10.5) 10"" , 



corresponding to 3.0 and 0.9 standard devia- 
tions, respectively. A graphical comparison of 
the results with the experimental value is given 
in Fig. 8. Also shown are our previous esti- 
mates [3,11] obtained before the CMD-2 and the 
new T data were available (see discussion below), 
and the recent evaluation of Hagiwara et al. [52] . 

10. DISCUSSION 

10.1. The Problem of the 27r Contribution 

The significant discrepancy between the e+e^ 
and T evaluations of ahad,LO jg ^ niatter of con- 
cern. The following changes in the dominant 
vr+TT^ contribution (all expressed in 10"^" units) 
are observed with respect to our earlier work [11]: 

• the new CMD-2 data [23] produce a down- 
ward shift of the e+e~ evaluation by 1.9 
(well within errors from previous experi- 
ments), 
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• the new ALEPH data [25] increases the r 
evaluation by 3.5, 

• including the CLEO data in the t evalua- 
tion slightly improves the precision, but fur- 
ther raises the central value by 4.0, 

• although including the OPAL data has lit- 
tle effect on the overall precision, it also in- 
creases the result by 1.9, 

• the new complete isospin symmetry- 
breaking correction, including the re- 
evaluation of the Sew factor, increases the 
r evaluation by 0.2. 

In principle, the observed discrepancy for the 27r 
contribution, (-21.2±7.3), or (-4.2±1.4)% when 
expressed with respect to e+e~, could be caused 
by any (or the combination of several) of the fol- 
lowing three effects which we examine in turn: 

• The normalization of e+e" data 
Here, as below, 'normalization' does not 
necessarily mean an overall factor, but 
refers to the absolute scale of the 'bare' 
cross section at each energy point. There 
is no cross check of this at the precision of 
the new CMD-2 analysis. The only test we 
can provide is to compute the e~^e~ inte- 
gral using the experiments separately. Be- 
cause of the limited energy range where 
the major experiments overlap, we choose 
to perform the integration in the range of 

from 610.5 to 820 MeV. The corre- 
sponding contributions are: 313.5 ± 3.1 for 
CMD-2, 321.8± 13.9 for OLYA, 320.8±12.6 
for CMD, and 323.9 ±2.1 for the isospin- 
corrected r data. No errors on radia- 
tive corrections and isospin breaking are in- 
cluded in the above results. 

• The normalization of r data 

The situation is quite similar, as the evalu- 
ation is dominated by the ALEPH data. It 
is also possible to compare the results pro- 
vided by each experiment separately, with 
the spectral functions normalized to the re- 
spective hadronic branching ratios. Leav- 
ing aside the region below 500 MeV where 



a fit combining analyticity constraints is 
used, the contributions are: 460.1 ± 4.4 
for ALEPH, 464.7 ± 9.3 for CLEO and 
464.2 ±8.1 for OPAL, where the common 
error on isospin breaking has been left out. 
The three values are consistent with each 
other and even the less precise values are 
not in good agreement with the e+e~ esti- 
mate in this range, 440.8±4.7, not including 
the error on missing radiative corrections. 

Apart from an overall normalization effect, 
differences could originate from the shape 
of the measured spectral functions. If all 
three spectral functions are normalized to 
the world average branching ratio (our final 
procedure), then the results for the contri- 
bution above 0.5 GeV become: 459.9 ± 3.6 
for ALEPH, 465.4 ± 5.1 for CLEO and 
464.5 ±5.1 for OPAL, with a common error 
of ±2.4 from the tttt^ and leptonic branch- 
ing ratios and the uncertainty on isospin 
breaking left out. Again the results are con- 
sistent and their respective experimental er- 
rors give a better feeling of the relative im- 
pact of the measurements. 

• The isospin-breciking correction ap- 
plied to T data 

The basic components entering SU(2) 
breaking have been identified. The weak 
points before were the poor knowledge of 
the long-distance radiative corrections and 
the quantitative effect of loops. Both points 
have been addressed by the analysis of 
Ref. [31] showing that the effects are small 
and covered by the errors previously ap- 
plied. The overall effect of the isospin- 
breaking corrections (including FSR) ap- 
plied to the 27r r data, expressed in relative 
terms, is (— 1.8±0.5)%. Its largest contribu- 
tion (—2.3%) stems from the uncontrover- 
sial short-distance electroweak correction. 
Additional contributions must be identified 
to bridge the observed difference. 

Thus we are unable at this point to identify 
the source of the discrepancy. More experimen- 
tal and theoretical work is needed. On the ex- 
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perimental side, additional data is available from 
CMD-2, but not yet published. As an alternative, 
a promising approach using e~^e~ annihilation 
events with initial state radiation (ISR), as pro- 
posed in Ref. [54], allows a single experiment to 
cover simultaneously a broad energy range. Two 
experimental programs are underway at Frascati 
with KLOE [55] and at SLAC with BABAR [56]. 
The expected statistics are abundant, but it will 
be a challenge to reduce the systematic uncer- 
tainty at the level necessary to probe the CMD-2 
results. As for r's, the attention is now focused 
on the forthcoming results from the B factories. 
Again, the quality of the analysis will be deter- 
mined by the capability to control systematics 
rather than the already sufHcient statistical ac- 
curacy. On the theory side, the computation of 
more precise and more complete radiative correc- 
tions both for e+e~ cross sections and r decays 
should be actively pursued. 

Other points of discussions are pursued in 
Ref. [22]. 

10.2. Consequences for a(M|) 

In spite of the fact that the present analysis 
was focused on the theoretical prediction for the 
muon magnetic anomaly, it is possible to draw 
some conclusions relevant to the evaluation of 
the hadronic vacuum polarization correction to 
the fine structure constant at M|. The problem 
found in the 27r spectral function is less important 
for Aa{M^) with respect to the total uncertainty, 
because the integral involved gives less weight to 
the low-energy region. The difference between the 
evaluations using the 27r, Aw and 27r27r'' spectral 
functions from e+e~ and r data are found to be: 

A<d(^l) - AaLd(^l) = (-2-8 ± 0.8) 10"^ . 

While this low-energy contribution shows a 3.5 
standard deviation discrepancy (when adding the 
different errors in quadrature), it also exceeds the 
total uncertainty of 1.6 10^^ on Aq:(71/|) which 
was quoted in Ref. [11]. It is worth pointing out 
that such a shift produces a noticeable effect for 
the determination of the Higgs boson mass Affj in 
the global electroweak fit [57]. With the present 
input for the electroweak observables [53] from 
LEP, SLC and FNAL yielding central values for 



Mh around 100 GeV, going from the e+e~ to the 
T-based evaluation induces a decrease of Mh by 
16 GeV using all observables and by 20 GeV when 
only the most sensitive observable, (sin^ Ow)eS, is 
used. 

11. CONCLUSIONS 

A new analysis of the lowest-order hadronic 
vacuum polarization contribution to the muon 
anomalous magnetic moment has been presented. 
It is based on the most recent high-precision ex- 
perimental data from e^e^ annihilation and r de- 
cays in the tttt channel. Special attention was 
given to the problem of isospin symmetry break- 
ing and the corresponding corrections to be ap- 
plied to T data. The main results of our analysis 
are the following: 

• the new evaluation based solely on e+e~ 
data is significantly lower than previous es- 
timates and is in conflict with the experi- 
mental determination of ttfj, by 3.0 standard 
deviations. 

• the new precise evaluations of the dominant 
TTTT contributions from e^e^ annihilation 
and isospin-breaking corrected r decays are 
not anymore in agreement with each other. 
A discussion has been presented for possi- 
ble sources of the discrepancy which could 
not be resolved. This situation is a matter 
of great concern, as the r-based prediction 
of is in better agreement with the ex- 
perimental value, from which it deviates by 
non-significant 0.9 standard deviations. 

More experimental and theoretical work is needed 
to lift the present uncertainty on whether or not 
new physics has been uncovered with the muon 
magnetic moment. 
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